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Abstract We present calculations of deuteron elastic and inelastic breakup cross sections and angular distri¬ 
butions at deuteron energies below 100 MeV obtained using the post-form DWBA approximation. The elastic 
breakup cross section was extensively studied in the past. Very few calculations of inelastic breakup have 
been performed, however. We also analyze the angular momentum - energy distributions of the cross section 
for formation of the compound nucleus after inelastic breakup. 
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1 Introduction 

Deuteron-induced reactions are being used to produce medical radioisotopes m and as surrogates to other 
reactions (see review El and references therein), among recent applications. Although they have been studied 
for decades ESI 13 [61, the complexity of these reactions continues to make their theoretical description 
challenging. 

The direct reaction mechanism is a major contributor to the deuteron reaction cross section due to the 
particle’s low binding energy. . Competition between elastic breakup, absorption of only a neutron or a proton 
(stripping and inelastic breakup) and absorption of the deuteron must be taken into account to determine the 
formation or not of a compound nucleus and its subsequent decay. The inelastic breakup reactions - those in 
which either only a neutron or a proton is absorbed - are particularly complex, as they form compound nuclei 
with a wide range of excitation energies and angular momenta. We present the results of a theoretical study of 
elastic and inelastic deuteron breakup for a large selection of targets at incident deuteron energies below 100 
MeV. We use the zero-range post-form DWBA approximation to calculate the elastic breakup cross section (3J 
@‘| and its extension to absorption channels to calculate the inelastic breakup cross sections EO- We discuss 
the regularities and ambiguities in our results, as well as the irregularities in the inelastic breakup energy and 
angular momentum distributions that complicate their substitution by a smooth distribution obtained from 
systematics. 
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2 Theory 

A reasonable theoretical description of elastic deuteron breakup was developed and applied to a multitude of 
experimental data almost forty years ago by G. Baur and collaborators (HH . The double-differential inelastic 
breakup cross section can be written in terms of its T-matrix element as 

(jbu 9 jr 1 

= 8 ( E d + £d-E p -E n ) ( 1 ) 

where k n and k p are the final neutron and proton momenta, respectively, kj is the initial deuteron momentum 
and the sum of neutron and proton final kinetic energies is constrained to the sum of the initial deuteron 
kinetic energy and its binding energy by the 8 -function. The T-matrix element can be well approximated by 
the post-form of the DWBA matrix element, 

r(k p ,k„;kj) = (\fr { p~\k. p ,r p )\jlrt\kn,rn) |v>(r)| v4 +) ( k rf> R )<fe( r )) , (2) 

which, in turn, can be well-approximated within the zero-range DWBA approximation 

T(k p ,k n -,k d )^D 0 (v^“ } (k p ,aR)v^ _) (k„,R)A(/?) | v4 +) (k d ,R)) (3) 

by including a correction for finite-range effects A (R) 0 and taking Dq =-125 Mev-fm 3 / 2 . 

To calculate the inelastic breakup cross section, one first analyzes the inclusive double differential breakup 
cross section, here given for the proton, again in the post-form of the DWBA. The initial state of the target is 
its ground state, <Pa, but the final neutron-target state, i \f c nA , can be any composite state allowed by energy and 
angular momentum conservation, 


d 5 o 2 % 1 


d)<? — fivd (2;r) 3 ^ \^p^ r p)VnA | v /w( r )| v4 S(Ed~\-£d E p E nA ). (4) 

Following Kasano and Ichimura 0 , we write the 5-function as the imaginary part of an energy denominator, 


d 3 a 

dk l 


hv d (2 7l) 3 


Im T(v4 +) ^^|v p „| y/j, } yfn A 

c 

x (E+ + £ d -E p -E c ^)~Uxjfj~ ) \if c nA \v pn \^ +) <l) d <P A 


the target ground-state matrix element of which we interpret as a neutron optical propagator, 

(+) ,„ - ' ^ / ^ ^ ^ 


G^(E d + £ d -E p ) = ^(<P A \V nA ) 


E d + £ d E p E c nA 


(VnA I &a) ■ 


(5) 


This furnishes a cross section of the form 

5ip-^(2^ In V" (r ' ) 

where the effective neutron wave function is given by 


d + \E d +£ d 



Xn( r») = (y { p \r p ) \v pn (r)\ v4 +) ( R )^(r)) • 

To reduce this further, the imaginary part of the optical propagator is decomposed as 

ImG„ = (1 + G\ul) ImG’o (1 + U n G n ) + G\W n G n . 


( 6 ) 


(7) 


( 8 ) 
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The inclusive proton emission cross section from breakup can then be separated into an elastic and inelas¬ 
tic part, denoted here by bu for elastic breakup and bf for inelastic breakup or breakup-fusion, corresponding 
to the first and second terms of the decompositon of the propagator, respectively, 


d 3 o _ d 3 o bu d 3 o bf 
dk 3 dk 3 + dk 3 


The contribution due to elastic breakup 


d 3 o bu 

dk l 


2K 1 f d 3 k n 

fadWfj J 2 n? 


T(k p , k n ;krf) | :2 8(E d + e d - E p 


E n ), 


(9) 


( 10 ) 


is just the double differential cross section of Eq. 0 integrated over the neutron momentum. The inelastic 
breakup cross section takes the form of an expectation value of the imaginary part of the optical potential, 


d 3 o bf 


2 1 
hv d (2 k) 3 


(kp,r n ,kj) | W n (v n ) | ^P n (k p: v n ;kj)^ , 


where the effective neutron wave function is given by 


(ID 


|fn(kp,r«;kj)) = \k p ,r p )G^ (r n ,r' n ) |v p „(r)| y4 +) (krf, R )<M r )) . (12) 

The physical interpretation of this cross section is simple: the deuteron first breaks up and, after propagating 
further, the neutron is absorbed while the proton is emitted. This wave function can be well-approximated in 
the zero-range approximation, again including the finite range correction A (R) of Ref. 0, as 


\%(k p ,r n ;k d ))^D 0 (y/j, } (k^aR)Gi +) (r n ,R)A(R) | ^ + \k d ,R)) . (13) 

To perform numerical calculations, the wave functions and matrix elements are expanded in partial waves 
of the orbital angular momentum alone. We thus neglect the effects of spin-orbit coupling but thereby reduce 
the number of matrix elements by a factor of almost 12. 

To calculate the distorted wave functions, the Koning-Delaroche global optical potentials [8) were used in 
the proton and neutron channels while the potentials of Refs. 0 or m were used to describe the deuteron 
scattering. Both deuteron optical potentials yield very similar results. We used the potential of Ref. (3 in the 
calculations presented here. 

The elastic breakup matrix elements of Eq. 0 are only conditionally convergent. A simple brute force 
application of standard integration methods requires integrating to radii of the order of a nanometer. The 
numerical integration can be reduced to radii of the order of picometers by using asymptotic expansions of 
the Coulomb wave functions to approximate the integral in the external region analytically. However, we 
have found the most efficient means of performing the integrals to be their extension to the complex plane as 
proposed by Vincent and Fortune OH- In this case, we can usually limit the numerical integration to several 
hundreds of fm and could probably limit it even more, were we to invest more effort in the evaluation of the 
Coulomb wave functions in the complex plane. 


3 Calculations 

We begin by presenting calculations of neutron and proton spectra from the elastic and inelastic breakup of 
30 MeV deuterons incident on 48 Ti and 181 Ta, shown in Fig. [I] In both cases, the inelastic breakup spectra are 
larger than the elastic ones. The elastic and inelastic spectra for neutron or proton both peak at about the same 
energy in both reactions with the proton spectra peaking at higher energy than the neutron ones. 

The difference in the the maxima of the neutron and proton spectra in Fig.[l]can be interpreted in terms of 
the Coulomb deceleration of the deuteron and posterior Coulomb acceleration of the outgoing proton. At the 
radius at which breakup occurs, R / m , we would expect the deuteron to have lost a kinetic energy of Ze 2 /Rb u , 
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Fig. 1 (a) Proton and neutron emission spectra from elastic and inelastic breakup spectra of 30 MeV deuterons incident on 48 Ti. 
(b) Proton and neutron emission spectra from elastic and inelastic breakup spectra of 30 MeV deuterons incident on 181 Ta. 


which is later recovered by the outgoing proton. On the average, we would thus expect to find for the outgoing 
neutron and proton kinetic energies 


E n 


Ze 

- £d 

Kim 

Ze 2 

-( Ed + -E d 

K bu 


(14) 


Interpreting the difference in the peaks as the energy difference at the most probable radius, we find the 
breakup radius Rt, u to be 10.6 fm for 48 Ti and 14.6 fm for 181 Ta. We compare these values to strong interaction 
radii of 6.0 fm for 48 Ti and 8.6 fm for 181 Ta, obtained by taking Ri nt « 1.25 (A 1 / 3 +2 1 / 3 ) fm, with A the target 
mass number. We conclude that the breakup occurs predominantly in the Coulomb field of the target. 

In Fig. [2j we compare our calculations with the experimental inclusive double differential proton cross 
sections for 56 MeV deuterons incident on 27 A1 and 118 Sn ca. The cross sections include protons from 
both elastic and inelastic breakup. The calculations agree well with the experimental data for the case of 27 A1 
although the calculations are shifted to slightly higher energies, probably due to the fact that their center-of- 
mass motion is not extracted correctly. This shift is much smaller in the case of 118 Sn. However, here the 
spectrum is underpredicted by about 20% at an angle of 9.5° and by a few percent at 13°. As the breakup at 
small scattering angles (corresponding to large impact parameters) is principally elastic, we suspect that more 
partial waves should be included in our elastic breakup calculation. 



Fig. 2 (a) Inclusive double differential proton spectra for 56 MeV deuterons incident on 27 Al. (b) Inclusive double differential 
proton spectra for 56 MeV deuterons incident on 118 Si. 
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Fig. 3 (a) Reaction probabilities as a function of the deuteron angular momentum for 25.5 MeV deuterons incident on 48 Ti. (b) 
Reaction probabilities as a function of the deuteron angular momentum for 80 MeV deuterons incident on 181 Ta. 


A standard optical model calculation of deuteron scattering furnishes a deuteron absorption probability 
(transmission coefficient) close to 1 below the grazing value of the angular momentum, corresponding to the 
solid black lines shown in Fig. [5] for deuterons incident at 25.5 and 80 MeV on 181 Ta. The flux lost due to 
elastic and inelastic breakup reduces this probability to about 75% for low partial waves in the nuclear interior 
and eliminates it in the surface region, as shown by the dashed lines in Fig. [3] Also shown are the neutron and 
proton absorption probabilities from inelastic breakup, which dominate in the surface region and extend into 
the nuclear interior with probabilities of about 10% and 20%, respectively. Only elastic breakup extends to 
higher angular momenta outside the range of the nuclear interaction. 

A complete calculation of a deuteron-induced reaction must also take into account the formation and 
decay of the compound nucleus (CN) formed by absorption of the deuteron, as well as those formed by 
aborption of the neutrons or protons of inelastic breakup. The deuteron absorption cross section is determined 
by the subtracted transmission coefficients, shown as dashed lines in Fig. [3] The excitation energy of the CN 
that results is determined by the kinetic energy of the incident deuteron. Inelastic breakup reactions produce 
protons and neutrons over the entire kinematically allowed range of energy. The corresponding differential 
formation cross sections are thus distributions in energy and angular momentum. Two examples of these, for 
neutron and proton absorption due to inelastic breakup in deuteron-induced reactions at incident energies of 
25.5 MeV and 80 MeV, respectively, are shown in Figs. [4] an d0 

In the reaction at 25.5 MeV, we observe in Fig. [4]that the neutron + target CN formation is well concen¬ 
trated at lower energies while the proton + target CN formation is concentrated for the most part at higher 
energies, due to the importance of the Coulomb repulsion at this energy. The separation in energy of the neu¬ 
tron and proton distributions is still visible but much less distinct at 80 MeV, as can be seen in Fig. [5] At both 



Fig. 4 (a) Energy - angular momentum distribution of neutron + target differential CN formation cross section for 25.5 MeV 
deuterons incident on 181 Ta. (b) Energy - angular momentum distribution of proton + target differential CN cross section forma¬ 
tion for 25.5 MeV deuterons incident on 181 Ta. 
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Fig. 5 (a) Energy - angular momentum distribution of neutron + target differential CN formation cross section for 80 MeV 
deuterons incident on 181 Ta. (b) Energy - angular momentum distribution of proton + target differential CN cross section forma¬ 
tion for 80 MeV deuterons incident on 181 Ta. 


incident energies, the proton + target distribution extends to larger values of the angular momentum, mainly 
due to the higher average energies of the protons that are absorbed. 

Well-defined structures in energy and angular momentum appear in the differential cross sections at the 
lower incident deuteron energy of 25 Mev and, to a lesser extent, in the distributions at 80 MeV. These 
reactions should thus be used with care as surrogates to other reactions 0. The compound nucleus formed 
could have an initial energy-angular momentum distribution very different from the one expected of a neutron- 
or proton-induced reaction. 


4 Summary 

We have used the post-form DWBA to calculate elastic and inelastic deuteron breakup cross sections in the 
zero-range limit. The breakup occurs for the most part outside the range of the nuclear interaction. However, 
the inelastic breakup cross sections, in which either the proton or neutron is absorbed by the target, are larger 
than the elastic one, in which the neutron and proton are simultaneously emitted. 

The breakup reactions reduce deuteron absorption at small impact parameters by about 25% and dominate 
the surface region completely. The neutron + target absorption of inelastic breakup tends to occcur at lower 
energy than the proton + target one, due to Coulomb repulsion, which inhibits low-energy protons from 
entering the range of the nuclear interaction. The low-energy inelastic breakup CN formation cross sections 
are found to have well-defined structures in energy and angular momentum, which could have important 
consequences in reactions in which they are intended to serve as surrogates. 
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